A model of driven and decaying magnetic turbulence in a cylinder 
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Using mean-field theory, we compute the evolution of the magnetic field in a cylinder with outer perfectly 
conducting boundaries, an imposed axial magnetic and electric field. The thus injected magnetic helicity in the 
system can be redistributed by magnetic helicity fluxes down the gradient of the local current helicity of the 
small-scale magnetic field. A weak reversal of the axial magnetic field is found to be a consequence of the 
magnetic helicity flux in the system. Such fluxes are known to alleviate so-called catastrophic quenching of the 
a effect in astrophysical applications. Application to the reversed field pinch in plasma confinement devices is 
discussed. 
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I. INTRODUCTION 

The interaction between a conducting medium moving at 
speed U through a magnetic field B is generally referred 
to as a dynamo effect. This effect plays important roles in 
astrophysics yj, |2[], magnetospheric physics [3], as well as 
laboratory plasma physics [4]. It modifies the electric field 
in the rest frame, so that Ohm's law takes the form J = 
a (E + U x B), where J is the current density, E is the elec- 
tric field, and a is the conductivity. Of particular interest for 
the present paper is the case where an external electric field 
E CKt is induced through a transformer with a time-varying 
magnetic field, as is the case in many plasma confinement ex- 
periments. With the external electric field included, Ohm's 
law becomes 

J = a (E + E c ^ + U x B) . (1) 

In a turbulent medium, often only averaged quantities (indi- 
cated below by overbars) are accessible. The averaged form 
of Ohm's law reads 

J = <t(e + if xt + U x B + 1\ , (2) 



where £ = u x b is referred to as the mean £ = u x b 
is referred to as the mean turbulent electromotive force, and 
u = U — U and b = B — B are fluctuations of velocity and 
magnetic field, respectively. It has been known for some time 
that the averaged profiles, J and oE do not agree in actual 
experiments. This disagreement cannot be explained by the 
U x B term either, leaving therefore £ as the only remain- 
ing term. Examples include the recent dynamo experiment in 
Cadarache JH] and in particular the reversed field pinch (RFP) 
0H01j which is one of the configurations studied in con- 
nection with fusion plasmas. The name of this device derives 
from the fact that the toroidal (or axial, in a cylindrical geom- 
etry) magnetic field reverses sign near the periphery. Indeed, 
in the astrophysical context it is well-known that the £ is re- 
sponsible for the amplification and maintenance of large-scale 
magnetic fields |[lj,|2|]. 



The analogy among the various examples of the £ term 
has motivated comparative research between astrophysics and 
plasma physics applications (8[]. In these cases, £ is found 
to have a component proportional to the mean field (aB, re- 
ferred to as the a effect) and a component proportional to the 
mean current density (r/tJ, where % is the turbulent diffusiv- 
ity). Since a is a pseudoscalar, one expects it to depend on 
the helicity of the flow, which is also a pseudoscalar. Deci- 
sive in developing the analogy between the a effects in astro- 
physics and laboratory plasma physics is the realization that 
a is caused not only by helicity in the flow (kinetic a effect), 
but also by that of the magnetic field itself [9] . This magnetic 
contribution to the a effect has received increased astrophysi- 
cal interest, because there are strong indications that such dy- 
namos saturate by building up small-scale helical fields that 
lead to a magnetic a effect which, in turn, counteracts the 
kinetic a effect This process can be described quan- 

titatively by taking magnetic helicity evolution into account, 
which leads to what is known as the dynamical quenching for- 
malism that goes back to early work of Kleeorin & Ruzmaikin 
lfl3ll . However, it is now also believed that such quenching 
would lead to a catastrophically low saturation field strength 
13, unless there are magnetic helicity fluxes out of the do- 
main that would limit the excessive build-up of small-scale 
helical fields I151 . This would reduce the magnetic a effect 
and thus allow the production of mean fields whose energy 
density is comparable to that of the kinetic energy of the tur- 
bulence llql . 

These recent developments are purely theoretical, so the 
hope is that more can be leamt by applying the recently gained 
knowledge to experiments like the RFP la, |7|] . Unlike toka- 
maks, the RFP is a relatively slender torus, so it makes sense 
to study its properties in a local model where one ignores cur- 
vature effects and considers a cylindrical piece of the torus. 
Along the axis of this cylinder there is a field-aligned current 
that makes the field helical. This field is susceptible to kink 
and tearing instabilities that lead to turbulence. It is generally 
believed that the resulting mean turbulent electromotive force 
£ is responsible for the field reversal lH ll7ll , The turbulence 
is also believed to help driving the system toward a minimum 
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energy state |[l"8ll . This state is nearly force-free and main- 
tained by E . This adds to the notion that the RFP must 
be sustained by some kind of dynamo process lfl9ll . In Carte- 
sian geometry such a slow-down has previously already been 
modeled using the dynamical quenching formalism 12011 . 

The RFP has been studied extensively using three- 
dimensional simulations fl9ll2ll423ll . which confirm the con- 
jecture of J. B. Taylor 111 811 that the system approaches a min- 
imum energy state. Additional understanding has been ob- 
tained using mean-field considerations ll24l l25ll . Both, here 
and in astrophysical dynamos there is an a effect that quan- 
tifies the correlation of the fluctuating parts of velocity and 
magnetic field. However, a major difference lies in the fact 
that in the RFP the a effect is caused by instabilities of the 
initially large-scale magnetic field while in the astrophysical 
case one is concerned with the problem of explaining the ori- 
gin of large-scale fields by the a effect JUS]- However, this 
distinction may be too simplistic and there is indeed evidence 
that in the RFP the a effect exists in close relation with a finite 
magnetic helicity flux ll26ll . supporting the idea that so-called 
catastrophic quenching is avoided by helicity transport. 

The purpose of this paper is to apply modern mean-field dy- 
namo theory with dynamical quenching to a cylindrical con- 
figuration to allow a more meaningful comparison between 
the a effect in astrophysics and the one occurring in RFP ex- 
periments. 



H. THE MODEL 

To model the evolution of the magnetic field in a cylinder 
with imposed axial magnetic and electric fields, we employ 
mean-field theory, where the evolution of the mean field B 
is governed by turbulent magnetic diffusivity and an a effect. 
Unlike the astrophysical case where a depends primarily on 
the kinetic helicity of the plasma, in turbulence from current- 
driven instabilities the a effect is likely to depend primarily 
on the current helicity of the small-scale field [9]. The current 
density is given by J = V x B j po, where po is the vacuum 
permeability and j = V x b/po is the fluctuating current 
density. The mean current helicity density of the small-scale 
field is then given by j ■ b. To a good approximation, the j ■ b 
term is proportional to the small-scale magnetic helicity den- 
sity, a ■ b, where a = A — A is the vector potential of the 
fluctuating field. The generation of a ■ b is coupled to the de- 
cay of A ■ B through the magnetic helicity evolution equation 
|[S Ell EH such that A ~B + a b evolves only resistively 
in the absence of magnetic helicity fluxes. 

Note that a ■ b is in general gauge-dependent and might 
therefore not be a physically meaningful quantity. However, if 
there is sufficient scale separation, the mean magnetic helicity 
density of the fluctuating field can be expressed in terms of 
the density of field line linkages, which does not involve the 
magnetic vector potential and is therefore gauge-independent 
ll28ll . For the large-scale field, on the other hand, the magnetic 
helicity density does remain in general gauge-dependent [29], 
but our final model will not be affected by this, because the 
magnetic helicity of the large-scale magnetic field does not 



enter in the mean-field model. 

We model an induced electric field by an externally applied 
electric field E cxt . In the absence of any other induction ef- 
fects this leads to a current density J — aE cxt . Furthermore, 
we ignore a mean flow (U = 0), and assume that the velocity 
field has only a turbulent component u. For simplicity we as- 
sume that E has no fluctuating part, i.e. E cxt = E . The 
decay of B is accelerated by turbulent magnetic diffusivity rjt, 
which is expected to occur as a result of the turbulence con- 
nected with kink and tearing instabilities inherent to the RFP. 
This mean turbulent electromotive force has two components 
corresponding to the a effect and turbulent diffusion with 

£ = aB - %p J, (3) 

where we have ignored the fact that a effect and turbulent 
diffusion are really tensors. The evolution equation for B is 
then given by the mean-field induction equation, 

^ = V x (aB - vthoJ + E * 1 ) , (4) 

where rfr = f] t + rj is the sum of turbulent and microscopic 
(Spitzer) magnetic diffusivities (not to be confused with the 
resistivity r/po, which is also often called 77). Note that only 
non-uniform and non-potential contributions to E can have 
an effect. 

As a starting point, we assume that the rms velocity it rms 
and the typical wavenumber fef of the turbulence are constant, 
although it is clear that these values should really depend on 
the level of the actual magnetic field. We estimate the value 
of r]t using a standard formula for isotropic turbulence, 

r) t = \tu 2 , (5) 

where r = (« rm s^f) 1 is the correlation time of the turbu- 
lence and M rms = (it 2 ) 1 / 2 is its rms velocity. Thus, we can 
also write rjt = u Tms /3k{. The turbulent velocity results from 
kink and tearing mode instabilities and will simply be treated 
as a constant in our model. For the a effect we assume that 
the kinetic helicity is negligible It would be much smaller than 
the current helicity, but of the same sign H, so it would con- 
tribute to quenching the a effect, and so we just take 

a = \rj ■ b/po, (6) 



and use the fact that j ■ b and a ■ b are proportional to each 
other. Here, po is the mean density of the plasma. For ho- 
mogeneous turbulence we have j ■ b = kfa ■ bj /xq, although 
for inhomogeneous turbulence, kfa-b/po has been found to 
be smaller than j ■ b by a factor of two (3(3]. We compute the 
evolution of a ■ b by considering first the evolution equation 
for A ■ B. Note that A ■ B evolves only resistively, unless 
there is material motion through the domain boundaries N29II . 
so we have 

n pvt 

— A B = 2E ■ B - 2r)p J • B - V • :F, (7) 
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where T is the mean magnetic helicity flux. While B evolves 
subject to the mean field equation the magnetic helicity of 
the mean field will change subject to the equation 

^-(A-B) = 2£ tot ■B-2r)u J-B- V -T m , (8) 

dt x ' 

where £ tot = £ + l° xt and F m = Ex is the mean 

magnetic helicity flux from the mean magnetic field, and $ 
is the mean electrostatic potential. Here, E = rj/ioJ — £ tot 
is the mean electric field. Subtracting (0 from (|7), we find a 
similar evolution equation for a ■ b, 

—a^b = -2£-B-2r)fJ, J 7 b- V-^ f , (9) 

where T{ = T — J- m is the mean magnetic helicity flux from 
the fluctuating magnetic field. Note that E cxt does not enter 

^ cxt 

in Eq. (O, because E cxt — E has no fluctuations. This 
equation can readily be formulated as an evolution equation 
for a by writing a = (rkf /3poA*o) a ■ b, i.e., 

= -2 Vt kf£ ■ B/B* q - 2 V kfa - V • F a , (10) 

where T a — (rkf /3poA*o which is a rescaled mag- 
netic helicity flux of the small-scale field and B cq is the field 
strength for which magnetic and kinetic energy densities are 
equal, i.e., 

#eq = MoPou? ms = (3poMo/t) Vt- (11) 

We recall that in the astrophysical context, equation (fTOt is 
referred to as the dynamical quenching model [fl~3l l27ll . In a 
first set of models we assume T a = 0, but later we shall 
allow for the fluxes to obey a Fickian diffusion law, 

T a = -/t Va, (12) 

where K a is a diffusion coefficient that is known to be compa- 
rable to or somewhat below the value of % 1291 [30ll . 

We solve the governing one-dimensional equations (|4]i 
and ([Tol l using the PENCIL CODE in cylindrical coordinates, 
(r, (f>, z), assuming axisymmetry and homogeneity along the z 
direction, d/dcf) = d/dz = 0, in a one-dimensional domain 
< r < R. On r = regularity of all functions is obeyed, 
while on r — R we assume perfect conductor boundary con- 
ditions, which implies that nxE = nxJ = 0, and thus 
n x dA/dt = 0, i.e., h x A = const. 

As initial condition, we choose a uniform magnetic field Bq 
in the z direction. In terms of the vector potential, this implies 

A(r,0) = (0,fl r/2,0) (13) 

for the initial value of A(r, t). 

We drive the system though the externally applied mean 
electromotive force in the z direction. We choose 

E c z Xt (r)=S c xt J (A :i r), (14) 

where E is the value of the mean electromotive force on 
the axis and k\R sa 2.4048256 is the rescaled cylindrical 



wavenumber for which E z (R) = 0, which corresponds to 
the first zero of the Bessel function of order zero, and thus 
satisfies the perfect conductor boundary condition on r = R. 
An important control parameter of our model is the non- 
dimensional ratio 

p-vf 

Q = E /vtkfBo, (15) 

which determines the degree of magnetic helicity injection. 
Other control parameters include the normalized strength of 
the imposed field, 

B = B /B cq (16) 

and the value of Lundquist number, 

£ = VA/rjkf, (17) 

which is a nondimensional measure of the inverse microscopic 
magnetic diffusivity, where va — -Bo/VMoPo is the Alfven 
speed associated with the imposed field. The Lundquist num- 
ber also characterizes the ratio of turbulent to microscopic 
magnetic diffusivity, i.e., 

ft = Vt/v = u rms /3r]kf = C/3B, (18) 

which we refer to as the magnetic Reynolds number. Note 
that, if we were to define the magnetic Reynolds number as 
Rm — u Tms /r)kf, as is often done, then 7Z = R m /3 would be 
three times smaller. Finally, the wavenumber of the energy- 
carrying turbulent eddies is expressed in terms of the dimen- 
sionless value of kfR. We treat kf here as an adjustable pa- 
rameter that characterizes the degree of scale separation, i.e., 
the ratio of the scale of the domain to the characteristic scale 
of the turbulence. In most of the cases we consider k{R = 10. 
In summary, our model is characterized by four parameters: 
Q, B, C, and kfR. In models with magnetic helicity flux we 
also have the parameter K a /vt- 

In addition to plotting the resulting profiles of magnetic 
field and current density, we also determine mean-field mag- 
netic energy and helicity, as well as mean-field current helic- 
ity, i.e., 

M m = (B 2 /2fi ), H m = (A-B), C m - (J-B), (19) 

where (.) = .rdr/(^R 2 ) denotes a volume average and 
the subscript m refers to mean-field quantities. Following sim- 
ilar practice of earlier work ifTTi bill , we characterize the so- 
lutions further by computing the effective wavenumber of the 
mean field, fc m , and the degree e m to which it is helical, via 

k m = fJ, C m /H m , e m = C m /2fc m Af m . (20) 

In the following we shall refer to e m as the relative mag- 
netic helicity. We recall that, even though A ■ B is gauge- 
dependent, for perfect conductor boundary conditions, the in- 
tegral J A- B dV is gauge-invariant, and so is then k m . Simi- 
lar definitions also apply to the fluctuating field, whose current 
helicity is given by 

Ct = (a)B e 2 q //i ryt- (21) 
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The magnetic helicity of the fluctuating field is then Hf = 
HoCf/kj. The magnetic energy of the fluctuating field can be ®® 
estimated under the assumption that the field is fully helical, 
i.e., (6 2 ) = k f \(a ■ 6)|, so that M f = \Q |/2fc f . We study both ' 60 

the steady state case where Q and B are non-vanishing, and 
the decaying case where Q = B = 0. In the latter case, we I 
monitor the decay rates of the magnetic field. 
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III. RESULTS o.O 0.1 0.2 0.3 0.4 0.5 



A. Driven field-aligned currents 

We begin by considering the case without magnetic helicity 
fluxes and take B = 1, C = 1000 (corresponding to TZ = 333) 
and fcf R = 10. The resulting values of fc m and e m are given 
in Table J] and the mean magnetic field profiles are compared 
in Fig. [T]for different values of Q. It turns out that, as we in- 
crease the value of Q, the magnetic helicity of the mean field 
increases, i.e. the product e m fc m increases, but the relative he- 
licity of the mean magnetic field decreases slightly, i.e., e m 
decreases. The value of fc m increases with Q, which means 
that the mean field will be confined to a progressively thinner 
core around the axis. Furthermore, the anti-correlation be- 
tween e m and k m is also found when varying B (see TableHII). 
jC, or TZ. This is demonstrated in Fig. [2] where we show that 
e m is in fact proportional to (fc m /fci) -1 / 4 and that the prod- 
uct en^fcm/fci) 1 ^ 4 is approximately constant, even though Q, 
B, or C are varied. This scaling is unexpected and there is 
currently no theoretical interpretation for this behavior. 

It is interesting to note that k m does not vary significantly 
with B, provided TZ is held fixed. However, for weak fields, 
e.g., for B = 0.1, the dynamics of the mean field is no longer 
controlled by magnetic helicity evolution, and the value of 
k m R has then dropped suddenly by nearly a factor of 2, and 
e m is in that case no longer anti-correlated with k m . This data 
point falls outside the plot range of Fig. [2] and is therefore not 
included. Also, if only L is held fixed, so that TZ varies with 
B, then fc m is no longer weakly varying with B, and varies 
more strongly in that case. 



TABLE I: Q dependence of k m and e m for B = 1, C = 1000, and 
k f R = 10. 



Q 



0.01 0.03 0.10 0.20 0.50 1.00 



k m R 2.76 3.44 4.63 5.26 6.49 7.20 
e m 0.95 0.91 0.84 0.82 0.78 0.73 
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FIG. 1 : Equilibrium profiles for three different driving strengths for 
B = 1, £ = 1000, and hR = 10. 



We must ask ourselves why the axial field component does 
not show a reversal in radius, as is the case in the RFP. Ex- 
perimental studies of the RFP provide direct evidence for a 
reversal. By comparing radial profiles of the axial current, 
Jii/c, with those of the axial electric field, En , one concludes 



that the mismatch between the two must come from the £ 



term 11171 132I1 . These studies show that E\\ < J\\/cr near the 
axis and En > J\\/cr away from it (assuming Bn > on the 
axis). Comparing with Eq. ©, it is therefore clear that £n 
must then be negative near the axis and positive near the outer 
rim. Turning now to dynamo theory, it should be emphasized 
that there are two contributions to £«, one from aB and one 
from —rjt^oJ; see Eq. (0. Let us therefore discuss in the fol- 
lowing the expected sign of £ n. Given that Q is positive, J B 
must also be positive, and therefore we expect a to be positive. 
If the mean magnetic field was really sustained by a dynamo, 
the a term would dominate over the r/ t term, but this is likely 
not the case here. Indeed, by manipulating Eq. (TTOb we see 
that, in the steady state without magnetic helicity fluxes, the 
equation for a takes the form 



TABLE II: B dependence of fc m and e m for Q = 0.1, TZ - 
k f R = 10. 



100, and 



B 



0.1 0.2 0.5 



1 



10 



k m R 1.80 3.33 3.50 3.99 3.42 3.31 3.25 
e m 0.51 0.91 0.89 0.84 0.89 0.89 0.89 



a 



KrhfMoJ ■ B/B, 



1 + TZB 2 / B"L 



(22) 



see, e.g., Ref. Illtl . However, as alluded to above, the relevant 
term entering £ is the combination Q; re( j = a—rj^poJ '■ B j B , 
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FIG. 2: (Color online) Dependence of e m on fc m for different sets 
of runs where either Q is varied (filled symbols), B is varied while 
keeping 1Z = 100 (red diamonds), B is varied while keeping L = 
100 (orange triangles), or £ is varied (blue squares). 




which is the reduced a. Inserting Eq. d22i > yields 



a rc d 



r] t ^ Q J ■ B/B 
1 + 11B 2 /B2 



(23) 



with a minus sign in front. The important point here is that 
a Tcd is indeed negative if J ■ B is positive. This means that 
we can only expect E\\ < J\\/<J, which is the situation in the 
RFP near the axis ll32ll . In order to reverse the ordering and to 
produce a reversal of the axial field, one would need to have 
an a effect that dominates over turbulent diffusion. Note also 
that for strong mean fields, a re d is of the order of the micro- 
scopic magnetic diffusivity. (This situation is well-known for 
nonlinear dynamos, because there a rc d and the microscopic 
diffusion term r]k m have to balance each other in a steady state 

Q.) 

We note in passing that fc m enters neither in Eq. (l22l nor 
in Eq. d23l ). However, k m does enter if there is a magnetic 
helicity flux and it affects the time-dependent case, as is also 
clear from Eq. ([Tot . Both cases will be considered below. 



B. Effect of magnetic helicity flux 

Next, we study cases where a diffusive magnetic helic- 
ity flux is included. In our model with perfectly conducting 
boundaries, the magnetic helicity flux vanishes on the bound- 
aries, so no magnetic helicity is exported from the domain, 
but the divergence of the flux is finite and can thus modify 



FIG. 3: Effect of magnetic helicity flux on equilibrium profiles of 
B z for Q = 0.1 (upper panel) and Q = 0.03 (lower panel). Note 
the field reversal at the outer rim in the latter case. 



the magnetic a effect. The same is true of periodic bound- 
aries, where no magnetic helicity is exported, but the flux di- 
vergence is finite and can alleviate catastrophic quenching in 
dynamos driven by the kinetic a effect B34I1 . 

In Fig. [3]we compare profiles of B z with and without mag- 
netic helicity flux. It turns out that the n a term has the effect 
of smoothing out the resulting profile of B z . More interest- 
ingly, it can lead to a reversal of B z at intermediate radii. For 
our reference run with Q = 0.1 (upper panel), the reversal 
is virtually absent at the rim of the cylinder. This is mainly 
because the pinch is so narrow; see Table UTTI However, when 
decreasing Q to 0.03, there is a clear reversal also at the outer 
rim (lower panel). However, decreasing Q to 0.01 does not 
increase the extent of the reversal. In none of these cases the 
field reversal is connected with a change of sign of a rc d- In- 
stead, a rc d is always found to be negative, even in the pres- 
ence of a magnetic helicity flux. Thus, the sign reversal of B z 
is therefore associated with a sign reversal of J z at the same 
radius. Nevertheless, the reversal is still not very strong with 
min(S z )/ max(B 2 ) —0.07, while in laboratory RFPs this 
ratio is typically -0.2 JH]. 
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C. Decay calculations 

Next, we consider the case of a decaying magnetic field 
in the absence of an external electric field. In that case all 
components of B must eventually decay to zero. The evolu- 
tion of the magnetic energy of the resulting mean and fluc- 
tuating fields is shown in Fig. [4] together with the evolu- 

tion of km and e m . At early times, when (B ) ^ B^ q , the 
energy of the large-scale magnetic field decays at the resis- 
tive rate A = —21]^. During that time, the energy of the 
small-scale field stays approximately constant: the magneti- 
cally generated a effect almost exactly balances turbulent dif- 
fusion and the magnetic field can only decay at the resistive 
rate. However, at later times, when (B ) <C B^ q , the en- 
ergy of the small-scale field decays with a negative growth 
rate A = —2r/kf, which then speeds up the decay of the en- 
ergy of the large-scale magnetic field to a rate that is about 
1.3 x iitkf n , where we have used the value k m R = 3.1 that is 
relevant for the late-time decay. This value is also that obey- 
ing Taylor's [18] postulated minimum energy state. Again, no 
reversal of the magnetic field is found, except in cases where 
there is an internal magnetic helicity flux in the system. 



IV. CONCLUSIONS 

The present work is an application of the dynamical 
quenching model of modern mean-field dynamo theory to 
magnetically driven and decaying turbulence in cylindrical ge- 
ometry. In the driven case, an external electric field is applied, 
which leads to magnetic helicity injection at large scales. Such 
a situation has not yet been considered in the framework of 
mean-field theory. It turns out that in such a case there is a 
weak anti-correlation between the actual value of magnetic 
helicity of the mean field, e m k m , and the relative magnetic 

helicity e m with e m ~ fcm 1 ^ 4 - This weak anti-correlation is 
found to be independent of whether Q, B, or C are varied. 
No theoretical interpretation of this behavior has yet been of- 
fered. In the decaying case, we find that the decay rate is close 
to the resistive value when the field is strong, i.e., B > B eq , 
and drops to the turbulent resistive value when the mean field 
becomes weaker. This behavior has also been found in earlier 
calculation of the decay of helical magnetic fields in Cartesian 
geometry ||2f3l . 

The original anticipation was that our model reproduces 
some features of the RFP that is studied in connection with 



TABLE III: Values of fc m and e m with and without magnetic helicity 
flux. 
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FIG. 4: (Color online) Evolution of (B )/i?cq, fem/fei, and e m for 
different values of B. Note that time is shifted by t,, which is the 
time when fc m attains its second maximum. In the top panel, the 
red lines indicate resistive decay rate of the large-scale field at early 
times, resistive decay rate of the small-scale field at late times, and 
the turbulent decay rate of the large-scale field at late times. 



fusion plasma confinement. It turns out that the expected field 
reversal that gives the RFP its name is found to require the 
presence of magnetic helicity fluxes. Without such fluxes, 
there is no reversal; see Fig. [3] Nevertheless, the reversal is 
rather weak compared with laboratory RFPs. This discrep- 
ancy can have several reasons. On the one hand, we have 
been working here with a model that has previously only been 
tested under simplifying circumstances in which there is tur- 
bulent dynamo action driven by kinetic helicity supply. It is 
therefore possible that the model has shortcomings that have 
not yet been fully understood. A related possibility is that the 
model is still basically valid, but our application to the RFP 
has been too crude. For example, the assumption of fixed val- 
ues of r)t and B cq is certainly quite simplistic. On the other 
hand, it is not clear that this simplification would really affect 
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the outcome of the model in any decisive way. A different 
possibility is that the application of an external electric field 
is not representative of the RFP. However, an important basic 
idea behind the present setup has been to establish a model 
as simple as possible, that could be tested by performing cor- 
responding three-dimensional simulations of a similar setup. 
This has not yet been attempted, but this would clearly consti- 
tute a natural next step to take. 
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